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A: NON-EQUILIBRIUM NCA 

The extension of the non-crossing approximation onto 
the Keldysh contour has been discussed in several pa- 
pers [ll, 01 • It is customary to neglect the bias voltage 
dependence of the conduction electron density of states 
(DOS) p(e), This is justified provided p is well approxi- 
mated by a constant in a region around the Fermi energy 
that is large compared to the applied bias voltage. When 
the DOS vanishes in a power-law fashion at or near the 
Fermi energy, this is no longer possible and the equa- 
tions have to be generalized appropriately. The full set 
of equations to be solved for the two-channel pseudogap 
problem becomes 
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for the pseudo-fermion. The DOS {pL and pn) of the 
two leads do not have to be identical. The bias voltage 
applied across the system is eV'^""' = pl — P-r, where ml 
and Mi? are the chemical potentials of the left and right 
leads. 




FIG. 1: The STM measurement of the magnetic adatom in 
graphene. The S = 1/2 magnetic adatom is located at the 
center of the honeycomb lattice of graphene. 



B: FANO-LINESHAPES 

An experiment reminiscent of the situation consid- 
ered by us has been performed recently, where magnetic 
adatoms on graphene where investigated via scanning 
tunneling microscopy (STM), see Ref. 13. Our analysis 
can be extended to include the current-voltage character- 
istics measured by an STM (see FigH]). In this case, one 
of the two fermionic leads represents the STM tip and it 
is necessary to explicitly allow for the different tunneling 
paths between the STM tip, the adatom and the sub- 
strate which will act as the second lead. An important 
difference between the STM setup and our analysis so far 
is that the STM tip is a good metal, e.g. a single-channel 
lead with constant DOS at its Fermi energy. We here 
will model it by a two-channel lead with constant DOS 
at its Fermi energy. This is justified provided the cou- 
pling between the STM tip and the system is small as 
the RG scaling equations for two and one-channel case 
are identical up to fourth order in the tunneling matrix 
element. 

The theory of STM on magnetic adatoms on a metal sur- 
faces has been worked out by Schiller and Hershfield Q 
and by O. Ujsaghy et al. [J]. The current is obtained 
from 
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where p^ip is the density of states of the STM tip and Pgfj- 
is an effective density of states probed by the STM and 
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FIG. 2: Fano-lineshapes with Fano parameter q = 10 for 
various values of F (in units of the half-bandwidth D). Here, 
we have set ed — —0.2D, /i = — O.ID. 



depends on two tunneling rates i/ and tc that parame- 
terize the hybridization strength of the STM tip with the 
magnetic adatom (tf) and the graphene leads {tc). The 
effective density of states p^g can be recast into 

Pcfj = -lm\tlGciR,R,e) + {td + tcVGc{R,Rad,€)) 



X Gadie)itd + tcV*GciRad,R,e)) 



(2) 



where V is the hybridization strength between the 
graphene electrons and the magnetic adatom, Gc (e) is the 
advanced local graphene electron Green function at the 
locus of the STM tip R and Gc{R, Rad, e)) is the advanced 
graphene electron Green function connecting the locus of 
the tip with the position of the adatom at Rad, and tc 
(td) is the tunneling matrix element between the STM tip 
and the substrate (magnetic adatom). Gadi^) is the ad- 
vanced Green function of the magnetic adatom that can 
be obtained from the pseudo-particle Green functions of 
section A. 

In the linear regime, the Fano lineshape is given by the 
differential conductance dl /dV\v^o, which turns out to 
be proportional to the effective density of states pef / (e) : 
dI/dV\v-,o'xpcff{e^V). 

dI/dV\v^o can be cast into the Fano lineshape where 
the Fano parameter q is given by 
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and can be treated as approximately constant in the en- 
ergy range of interest [5] . 

Typical Fano-lineshapes in the linear regime are shown 
in Figjp] The 2CK behavior seen in the STM measure- 
ment j7| for Co-adatom at the center of the honeycomb 
lattice is signaled by the Kondo peaks at a; = /i in Fano- 




FIG. 3: Non- linear conductance for a large F — O.ID between 
impurity and graphene substrate and a much smaller hoping 
Tts between tip and the impurity, Vst <C F, corresponding to 
the STM measurement reported in [3. The G{V,T) curves 
agree well with the STM results of [31 • (a) logV dependence 
around V ~ T^. (b) behavior for V <T. (c) T^/^ 2CK 
behavior for T < V < Tk- Here, T = 5 x IQ-'^D. The other 
parameters are: p — —O.ID, — —0.2D. 



lineshapes, which are compatible with a large fitting pa- 
rameter q (for example q = 10) and a correspondingly 
small tc/td and concomitantly small intervalley scatter- 



ing. 



C. UNIVERSAL 2CK-LM CROSSOVER FOR p>0 

In the main text, we focus on the universal 2CK-LM 
crossover for negative chemical potential, /i < 0. A sim- 
ilar scaling behaviors can also be found in conductance 
for positive p. As shown in Fig. [H for a fixed positive 
p = 0.1-D, the linear conductance G{T) vs. hybridization 
r follows a single universal scaling form of T/T* . The 
single scaling form of G{T) we observe here for /i > is 
somewhat surprising as for ^ < the conductance shows 
two distinct scaling regimes: T < Tk and T > T*. We 
believe that this difference maybe due to the particle-hole 
asymmetry in our model as the Kondo peak, located at 
id = p, is affected more by the charge peak at < 0. for 
p < than that for p > 0. 

Similar to the case for p < 0, the linear conductance for 
p > shows a typical 2CK y/T behavior for T < Tk, and 
a universal power-law behavior at high temperatures for 
r -> r*: G(T) oc {T/T*Y with a « 0.00009 « 0. The 
Kondo temperature Tk and the crossover scale T* for 
/i > behave in a similar way to their p < Q counter- 
parts: Tk oc Fxe-i/r, T* cx {T-T*f/'' withF* w 0.061? 
and V ~ 0.05. We believe that our results for both posi- 
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FIG. 4: Universal scaling in linear conductance G{T) as a 
function of temperature T for at fixed positive chemical po- 
tential fj. — O.ID and for various values of hybridization F. (a) 
T^/^ 2CK behavior for T < Tk for large values of V. Inset: 
Tk/V vs. 1/r. (b) Universal scaling of G{T) for smaller val- 
ues of r. Inset: The crossover energy scale T* as a power-law 
function of F — F*. Here, ed = — 0.2_D, Go is the linear con- 
ductance for F = 0.62D at T = 5 x and F* « 0.06-D. 



tive and negative values of ^ could be used as theoretical 
guidance in future experiments to clarify the issue on 
two-channel Kondo physics in graphene. 
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Based on the non-crossing approximation, we calculate both the linear and nonlinear conductance 
within the two-lead two-channel single-impurity Anderson model where the conduction electron 
density of states vanishes in a power-law fashion oc |aj — /ifI' with r = 1 near the Fermi energy, 
appropriate for an hexagonal system. For given gate voltage, we address the universal crossover from 
a two-channel Kondo phase, argued to occur in doped graphene, to an unscreened local moment 
phase. We extract universal scaling functions in conductance governing charge transfer through the 
two-channel pseudogap Kondo impurity and discuss our results in the context of a recent scanning 
tunneling spectroscopy experiment on Co-doped graphene. 

PACS numbers: 



Introduction. The two-channel Kondo [ll (2CK) prob- 
lem is a fascinating example of an exotic quantum 
many-body phenomenon resulting in a metallic ground 
state with non-Fcrmi-liquid behavior. It involves a single 
quantum impurity spin with .s — 1/2 that couples antifer- 
romagnetically to two identical conduction electron reser- 
voirs. As a result, Kondo processes involving both reser- 
voirs lead to overscreening of the local moment. Theo- 
retically, the 2CK physics has been studied extensively 
via Bethe ansatz Q , conformal field theory Q , bosoniza- 
tion Q and the numerical renormalization group Q . Ex- 
perimentally, however, up to date, only very few exam- 
ples of clear 2CK physics have been experimentally real- 
ized, e.g. in semiconductor quantum dots in magnet- 
ically doped nanowires, and in metallic glasses [1, [J. 

Recently, the Kondo effect of magnetic adatoms in 
graphene has attracted much attention, theoretically [lol . 
as well as experimentally 



as well as experimentally [1^ due to the possible re- 
alization of a 2CK ground state. One interesting aspect 
of Kondo physics in graphene is due to the Dirac (linear) 
spectrum that gives raise to a pseudogap local density of 
states (DOS), Pc('j-') oc |a;|'' with r = 1, at the impurity 
site, making graphene one of the few experimental real- 
izations of the pseudogap Kondo model [ij, , the sim- 



(a)^ 



plest model to study critical Kondo destruction 1^ 16| . 
In the pseudogap Kondo (or more generally Anderson) 
model, a quantum phase transition is expected between 
the Kondo screened and the unscreened local moment 
(LM) ground states for < r < 1 [ll, [T 
the case of graphene, Kondo screening 
0," 



cur 
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For r = 1, 
does not oc- 
resulting in a LM ground state. Yet, 
Kondo screening can be induced by changing the Fermi 
energy, e.g. by applying a gate voltage (^ 7^ 0). 

That two independent screening channels can exists in 
graphene is related to the existence of two inequivalent 
Dirac points {K and K') in its band structure. It there- 
fore was argued that the effective low-energy model for 



LM 



T'~\'V 



universal 
crossover scaling 



2CK 



Tk 



iCK 
scaling 




FIG. 1: (a) Schematic phase diagram for the 2CK-LM 
crossover. The parameter j refers to either F or /i (in units of 
D — 1), and j* refers to the crossover scale for a fixed tem- 
perature To ~ 5 X 10~^-D. Tk and T* represent crossover en- 
ergy scales associated with the universal scaling for 2CK (blue 
shaded) and the high-temperature 2CK-LM crossover scaling 
(green shaded) regime, respectively, (b) Schematic setup of 
the model in Eq. (1): a S = 1/2 impurity (red dot) couples 
symmetrically to the two sub-lattices of two graphene leads 
that are kept at different chemical potentials of the setup, (c) 
Dirac points (labeled by K and K') with its Dirac spectrum 
on graphene's honeycomb lattice in momentum space. 



magnetic impurities in graphene depends on the loca- 
tion of the adatom [2l|: if the impurity is located at 
the center of the cell, the inter-valley scattering docs 
not couple the two screening channels and an effective 
2CK ensues [13 ■ This has been explicitly demonstrated 
within a tight-binding description where the hybridiza- 
tion between electronic states in graphene and irnpu- 
rity states preserves the A-B sub-lattice symmetry j22| . 
The situation is different if the adatom is located on a 
graphene site and the sub-lattice symmetry is absent in 
the impurity-graphene hybridization. The effective low- 
energy model is in this case the more conventional one- 
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FIG. 2: The impurity spectral function p(a;) vs. chemical 
potential fi (in units of the bandwidth D). The Kondo peak 
is pinned near jj,. The Dirac point is assigned oj = 0. Inset; 
Kondo temperature Tk (defined in text) as a function of ^. 



The NCA parameter are T 
-0.2D. 



5 X 10-^ D, r = 0.2D, ed = 



channel pseudogap Kondo model |2l|, |2^- The differ- 
ence between the various adatom positions can be p robed 
by scanning tunneling microscopy (STM) [23, 24|. This 
has been achieved recently with Co-atoms as magnetic 
adatoms where signatures of 2CK physics are seen when 
the adatom is located at the zone center of the honey- 
comb cell such that the inter-valley scatterings between 
K and K' Dirac electrons are strongly suppressed 12|. 

The 2CK-LM quantum phase transition in the two- 
channel pseudogap Kondo and Anderson model with 
fi = and < r < 1 in and out of equilibrium has 
been studied recently [li, [UHl, but the experimen- 
tally more relevant situation of r = 1 and fi ^ has 
not yet been properly addressed. The possibility of re- 
alizing 2CK physics was first pointed out in Ref. [lol |. 
There is, however, a lack of systematic investigations 
beyond the mean-field treatment. In this letter, we 
study the crossover phenomenon from the LM to 2CK 
regime, and work out the universal scaling functions of 
the conductance. A comprehensive analysis of the non- 
equilibrium transport, including STM lineshapes in the 
various regimes is presented. We address this issue using 
the non-crossing approximation (NCA) 2b 2^ : which is 
known to give reliable results for multi-channel Kondo 
systems which are in line with conformal field theory 
results 0. It has been shown recently that the NCA 
is able to reproduce the correct qualitative features of 
the two-channel pseudogap model for < r < 1 in and 
out of equilibrium in excellent agreement with exact re- 
sults 

Here, we study the steady-state transport properties of 
the experimentally most relevant pseudogap case, r = 1 
as a function of doping, temperature and bias voltage. 

The Model Hamiltonian. Our starting point is the 
two-lead two-channel single-impurity pseudogap Ander- 



son model. Each lead is characterized by a power law 
density of states (DOS) (pc(w) ~ |a;|9(£' — |a;|), where D 
is a high-energy cutoff that serves as our unit of energy 
D = 1). A possible realization of this model is sketched 
in Fig. 1 (b) where a spin- 1/2 impurity is coupled to two 
graphene leads with couplings to each lead preserving the 
sub- lattice symmetry of the honeycomb lattice 12|, |22| . 

Within a pseudofermion representation, the Hamilto- 
nian reads KTI 



k.a,T,a (7 



(1) 



where is the chemical potential of lead a, and a — 
L/R labels the two conduction electron leads. The in- 
dices (7 and T refer to spin and charge (related to K 
and K ) channels characterizing the conduction electrons. 
The second and third term of the RHS of Eq. ([T]) rep- 
resent the spin- 1/2 and the hybridization strength Ua 
between the graphene electrons and the impurity. In the 
pseudofermion representation the local electron is decom- 
posed as dl ^ ~ fa^T- Eq. ([T|) is a faithful representa- 
tion of the two-channel single-impurity pseudogap An- 
derson model, provided the constraint Q — fa-fo- + 
J2t ^t^t = 1 is fulfilled at all times. A finite bias voltage 
is implemented by shifting the chemical potentials in the 
leads such that ijll — I^^r = eV is the applied bias voltage 
across the 2CK system Q ■ 

To study the properties of Eq. ([J) we employ the NCA. 
Its ability to correctly capture the properties of the pseu- 
dogap two-channel Anderson model was established re- 
cently miii. Within the NCA, the retarded self-energy 
for pseudofermions, G^{uj) — [w — — S''(w]l~\ and 
slave-bosons, ^'■(w) = -W {uj)]-^ , are [ll[ll[2i,[27| 



E''(w) = -Y,j deT^{u - t - M„)/(e - uj - fi^,) {e) ,{2) 

The NCA expressions for the lesser self-energy of the 
pseudofermion, G'^{uj) = I]<(cj)|G"'(a;)p, and slave- 
boson, D<{u}) = Il<{uj)\D''{uj)\'^, are 

S<(w) = - E / deT^iiu - e - ti^)f{Lu ~ e ~ n^)D<ie),iA) 

II<{lu) = - E / deV^ie - Lu - iJ^)f{Lu - e + ti^)G<{e).i5) 

Here, Ta{uj) = ^aPc,a{^) with Fq, = ttIJ/q.^ with 
Pc.a{ijj) oc — fial SLud f (oj) = [1 + e^"]^"'" is the Fermi 
function. 
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The physical spectral function, p{uj,V), is the convo- 
lution of pseudofermion and slave-boson Greens function 



de[ImD''{e)G<{LO + e)- 

D<{e)ImG''{LJ + e)]. (6) 



The normalization factor Z ~ ^ J duj[D'^{uj) 
— G^(cj)] enforces the constraint, < Q >= 1. The cur- 
rent is given by 



I[V,T) = — / d(^—————p{uj,V,T) 
h J Tl(lu) + ThIlo) 

X [f{u + eV/2)^ f{uj-eV/2)]. 



(7) 



The nonlinear conductance G{V) = dI{V)/dV is com- 
puted by numerical derivative of the current I{V), 
whereas the linear-response conductance is directly ob- 
tained from 



G(0,T) 



h 

X p{uj,V 



duj 



2Tl{u)Tr{u:) 
0). 



df{Lo) 
dijj 
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Eqs.Q-® together with the Dyson equation for G(aj) 
and D{uj) forms a self-consistent set of integral equations. 
These equations are iterated until a solution is found with 
which Eqs.(l6l)-(l8]) can be evaluated. 

Results. Wc now turn to a discussion of the self- 
consistent solution of Eqs. (HI -([5]); the results are sum- 
marized in the phase diagram Fig. [Ija). For simplicity, 
we focus on the particularly simple case with parity (left- 
right) symmetry, Ul = Ur, Tl = Tr = ^. 

For larger values of j — T, \fi\ (or j — j* > O.ID) with 
r* ~ 0.05D and fi* ^ — 0.05D being crossover scales at 
a fixed temperature Tq ^ 10~^Z?, our results show clear 
2CK behaviors at low temperatures T < Tk where Tk is 
the Kondo temperature defined as the temperature where 
G(0, T) deviates from a ~ VT behavior |27j. in agreement 
with its conventional definition: G(0,Tr-) = G(0, 0)/2. 

However, for smaller values of j with j ~ j* < O.ID, 
we find universal power-law scaling of G{T, 0) distinct 
from both 2CK {i.e. ~ VT) and one-channel Kondo 
{i.e. ~ T^) [Hi behavior at temperatures T :$> T* where 
T* is the crossover energy scale, describing the 2CK-LM 
crossover. Note that the NCA gives reliable results for 
I{V, T) even in the single-channel Kondo case as long as 
T > O.ITk. 

The crossover scale is finite for any non- vanishing gate- 
voltage {fi 7^ 0). We checked that the crossover scales 
F*, vanish as T — in a power-law fashion, consis- 
tent with the general expectation that Kondo screening 
in graphene can be induced by arbitrarily small doping 
{p^O). 

The local density of states V), given by Eq.®, is 
shown in Figl2l The Kondo peaks occur at the chem- 
ical potential /i, and the Kondo temperature Tk fol- 
lows the pseudogap Kondo behavior with r = 1: Tk ^ 




(eV/k,T,)- 



FIG. 3: (a) (Go - G(0,T))/Go calculated from Eq.® dis- 
plays T^/^ behavior for T < Tk- Go = G(0, ro)r=o.225D is the 
linear conductance at the lowest numerically accessible tem- 
perature To « 5 X 10~^D at r = 0.225D. Inset: Exponential 
behavior for Tk ~ De"'"''''!'"". (b) An additional power law 
behavior at high temperature for small F « F* « 0.05D, 
where the crossover scale T* shows a power law behavior 
T* ~ IF - F*!^/" with u ~ 0.1 (Inset). Here, ea = -0.2D, 
ji — —O.ID. (c) ~ (e): Non-linear conductance at a large 
F = 0.2D, and fixed parameters: To = 5xlO"''D, = -O.ID, 
ed = -0.2D. The G(V, T) curve shows 2CK Kondo behavior: 
(c) logV dependence around V ^ Tk with Tk ~ 5 x 10~^D 
and T = To. (d) V'^ behavior for V < T. (e) T^/^ gCK be- 
havior for T < V < Tk. Here, the red lines are fits to the 
corresponding power-law or logarithmic behaviors in different 
bias regime. 



1^1 X e*^~°/l^l^ (Inset in Figl2]) where a has the unit of en- 
ergy and is a function of the bandwidth D and the Kondo 
couphng The \u}\ decay in the vicinity of u; = 

is a reflection of the Dirac point of the conduction elec- 
trons with DOS /9c(w) ~ jwl. Comparable results have 



also been discussed in Refs. 

Figs. [IKa) and [SKb) shows the linear response con- 
ductance G(0, T) obtained via numerical derivative of 
Eqs. ([7]) and ([8]). As long as ^ 7^ 0, we find clear 
2CK behavior at low temperature regime (Fig. [3ja)), 
with G(0, T) displaying a T^/"^ behavior for each coupling 
F [231 ■ The Kondo temperature Tk, behaves as: Tk ^ 
£)g-T|ed|/r (Fig, [3](a) Inset). For higher temperatures 
(Fig. [2Ib)), we find a universal power-law behavior in 
T/T* for T >T* near F* « 0.051): G(0,r) oc {T/T*^ 
with a « 0.029, indicating the universal crossover from 
2CK to LM regime. The crossover energy scale T* shows 
a power-law dependence, T* ~ |F - F*^/'' (Fig. [Sjb) 
Inset), where v ^ 0.1. 

G(y, T) at larger F shows clear 2CK behavior (see 
Fig. He)- (e) for F = 0.2D). Fig. Eljc) indicates the logV 
behavior predicted for Kondo scattering processes. The 
y2 behavior aXV <T and l/^/^ behavior aXT <V <Tk 
in Fig. [3]^d) and (e), respectively, are the characteristics 
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FIG. 4: Temperature scaling of the nonlinear conductance 
G{V, T) in the universal crossover region with F = 0.071) 
and = — 0.08-D (e^ = — 0.2D). All curves collapse onto a 
universal scaling function. It displays V^^'^ behavior at T < 
V < lOr (Inset). For large V (in units of D), V/T > 10^ 
G{V, T) shows a power law behavior similar to Fig[3jd). 




1 2 3 
(eV/kJj' 

FIG. 5: (a) Scaling of GiV, T) for large F (parameters are 
the same as in Fig[2]). For V < Tk, G{V,T) collapses onto a 
single curve ~ T^/^ (b) Power law behavior of G{V,T) for 
V > V* . Here, G{0, To, Fq) refers to the linear conductance at 
a fixed temperature To = 5 x 10~^i? and a fixed coupling Fq = 
0.2D. Inset: Power law behavior of the crossover energy scale 
V* r^\r - T'l" (see text). Here, ea = -0.2D, = -O.ID. 



of 2CK physics. 

To illustrate the crossover between the 2CK and LM 
regimes, G{V,T) is shown in Fig. |4]and Fig. [3 Fig. IH 
demonstrates V/T scaling [H of G(y,T) for F = 0.07L» 
and fj, = — 0.08D. The resulting curves collapse onto 
a universal function. For V/T > 10^, G{V,T) shows a 
power law behavior similar to the one shown in Fig. [HJb) ; 
for 1 < V/T < 10, the behavior follows the 2CK V^^^ be- 
havior as shown in the inset of Fig.|4l for V -^T/it shows 
V'^ Fermi-liquid behavior. Fig. [Sfa) is the 2CK scaling 
plot for large F, where all curves follow the scaling func- 
tion [13, G(0,T) - G{V,T) = By:T^/^H{A^). Here, 
H{x) oc a;^/^ is a universal function and and A are 
non-universal constants 0. For V > V*, G{V,T) shows 
the same universal 2CK-LM power-law crossover behav- 
ior for small coupling F w F* as shown in the linear con- 
ductance, G'(0,T): G{0,T) - G{V,T) oc {V/V*)°' with 
a w 0.029, (see Fig. [5{b)). The crossover energy scale 
V* also displays a power law behavior V* \T — F*]^/", 
where ~ 0.1, in line with the equilibrium behavior. 

Discussions and Conclusions. To make contact with a 
recent STM experiment [l^ , it is necessary to generalize 
our setup to the case where one of the leads is made up of 
a simple (i.e. one-channel) metal with constant DOS pup 
near its Fermi energy. The ground state in this case will 
be that of a Fermi liquid. If, however, the tip is coupled 
only very weakly (Ftip/Z) <C 1) to the magnetic adatom 
on graphene, the corresponding energy scale will be van- 
ishingly small. In this case it is permissible to replace the 
normal metal lead by a two-channel lead (with constant 
DOS at the Fermi energy), as the RG equations for the 
one-channel and two-channel Kondo problem coincide in 
lowest order. The behavior of G{V, T) in this case are 
compatible with the results discussed above, see Figs. 1 
and 2 of Ref. [2J|. Our results, see e.g. Fig. [3{c)-(e), are 



thus in line with the experimental findings reported in 
Ref. 12|. Fano-lineshapes, derived from our results are 



shown in 24 1. 

Universal out-of-equilibrium scaling is currently pursed 
in a wide range of condensed matter systems. As demon- 
strated, magnetic adatoms on graphene offer the pos- 
sibility to study steady state properties in a universal 
crossover regime. Here, we have addressed the univer- 
sal crossover regime that separates the local moment 
regime from the two-channel Kondo regime of adatoms 
in graphene. Our analysis is based on the two-channel 
Anderson model in graphene. In particular, we calcu- 
lated the differential conductance both in the linear and 
nonlinear regime. For sufficiently large hybridization, we 
found clear two-channel Kondo signatures. As the hy- 
bridization is reduced, the crossover region separating 
two-channel Kondo and local moment ground states is 
entered and the crossover is monitored by the narrowing 
of the Kondo resonances. In the crossover regime, the 
conductance shows universal power-law behavior. Our 
results not only agree well a recent scanning tunneling 
spectroscopy experiment but also provide a comprehen- 
sive theoretical analysis of the transport properties of 
2CK impurities in graphene. 
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